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In this paper, we study transport properties of non-equilibrium systems under the appHcation of 
hght in many-terminal measurements, using the Floquet picture. We propose and demonstrate that 
the quantum transport properties can be controlled in materials such as graphene and topological 
insulators, via the application of light. Remarkably, under the application of off-resonant light, 
topological transport properties can be induced; these systems exhibits quantum Hall effects in the 
absence of a magnetic field with a near quantization of the Hall conductance, realizing so-called 
quantum Hall systems without Landau levels first proposed by Haldane. 



I. INTRODUCTION 

Application of light is a powerful method to change ma- 
terial properties. For example, light can induce currents 
through mechanisms such as photovoltaic effect^, photo- 
thermoelectric effect^, and photo-drag effects^. More- 
over, light can change the response of materials and in- 
duce insulator to metal transitions^ or change the char- 
acteristics oi p — n junctions^! 

In recent years, there has been tremendous devel- 
opments and interests in the induction of quantum 
phases through light applications. For example, exper- 
iments have demonstrated that superconductivity can 
be induced through infrared pulses in high-temperature 
cuprate superconductors^. Inductions of quantum phases 
are inherently non-equilibrium phenomena, and thus 
their understanding is quite challenging. Even some ba- 
sic questions such as the physical signatures of the in- 
duced phases and how such phases can be stabilized in 
a steady state do not have answers yet. Many of quan- 
tum phases manifest themselves through transport, and 
therefore, understanding of transport properties in non- 
equilibrium, open systems is crucial for experimental ver- 
ifications of such induction of quantum phases under the 
application of light. In this paper, we develop a general 
formalism for studying non-equilibrium transport under 
the application of light, and, using the formalism, address 
the possibilities of the induction of topological properties 
through light. 

Motivated by recent rapid development of the under- 
standing in topological phases, the possibility of inducing 
topological phases such as integer quantum Hall phase 
and topological insulators through light has been theo- 
retically explored by many different groupS^HII'. Gener- 
ally speaking, the application of light on electron sys- 
tems has two important physical effects; 1. photon- 
dressing of band structures through the mixing of dif- 
ferent bands 2. redistribution of electron occupation 
numbers through the absorptions/emissions of photons 
leading to non-equilibrium distributions. Previous works 




FIG. 1: The many-terminal measurements of DC current for 
graphene under the application of light. Graphene is attached 
to multiple leads labeled by {a}, and the leads are connected 
to reservoirs at chemical potentials, {/i^}. Off-resonant, circu- 
larly polarized light is applied to the graphene. In the absence 
of impurities and interactions, electrons coming from leads co- 
herently propagate under the application of light, asborb or 
emit photons, and leak out into leads. Current measurements 
between each leads determine longitudinal and Hall conduc- 
tances. 



proposed optical induction of band structures with topo- 
logical propertie^fl^, and thus have mostly focused on 
the analysis of the first effect. On the other hand, most of 
these works do not address the question of the second ef- 
fect, the redistribution of electrons in the band structure, 
and thus its physics is yet poorly understood. Topolog- 
ical properties only appear when certain bands are fully 
filled, and it is not clear how this band occupation can 
be achieved and topological properties survive when the 
system is strongly driven out of equilibrium by the ap- 
plication of light. 

In order to answer these questions, we study the phys- 
ical consequence of the application of light through DC 
many-terminal transport measurements as in Fig.[lj The 
coupling of the driven systems with leads, that are in re- 
turn coupled with equilibrium reservoirs, plays the cru- 
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cial role to determine the occupations of electrons. Us- 
ing the formalism for transport properties in pe riodically 
driven systems developed by various groups^, here we 
study topological transport phenomena in materials such 
as graphene and three dimensional topological insula- 
tors. First of all, we show that non-equilibrium transport 
properties cannot generally be captured by the photon- 
dressed, effective band structures. In particular, in addi- 
tion to the usual transport through such static effective 
band structures, there are contributions from photon as- 
sisted electron conductions. Thus, the induction of topo- 
logically non-trivial hand structures does not immediately 
imply the topological properties of the non- equilibrium 
systems. 

On the other hand, the regime exists in which topolog- 
ical band structures can manifest themselves; we explic- 
itly demonstrate that for off-resonant light where elec- 
trons cannot directly absorb photons, the transport prop- 
erties of the non-equilibrium systems attached to the 
leads are well approximated by the transport properties 
of the system described by the static effective Hamil- 
tonian that incorporates the virtual photon absorption 
processes. In particular, the occupations of the electrons 
under this situation are close to the filling of the photon- 
dressed bands. As examples, we show that the transport 
properties under the application of off-resonant light is 
given by the photon-dressed Hamiltonian corresponding 
to a quantum Hall insulator without Landau levels^ in 
the case of graphene, and to a gapped insulator with 
anomalous quantum Hall effects and magneto-electric re- 
sponse described by axion electrodynamics^ at the sur- 
faces of three dimensional topological insulators. In these 
systems, the measurements in six-terminal configurations 
in Fig. [l] lead to the near quantization of Hall conduc- 
tance. Thus, the application of off-resonant circularly 
polarized light leads to an intriguing "Hall" effects with- 
out applying a static magnetic field. 

This paper is organized as follows. In Section [Hj we 
describe the summary of the results, focusing on the anal- 
ysis of graphene and three dimensional topological insu- 
lators under the application of off-resonant light. Here 
we provide the physical and intuitive explanations of the 
phenomenon of light-induced quantum Hall effects and 
refer to later sections for many important details. In Sec- 
tion [TTTj we develop the formalism for studying the non- 
equilibrium transport properties under periodical drives 
in many-terminal measurements. Our formalism is based 
on the extension of the multi-probe Biittiker-Landauer 
formula^^ to periodically driven systems, a "Floquet Lan- 
dauer formula." We provide two distinct ways to calcu- 
late the transmission amplitudes in the driven systems. 
First method expresses the results in terms of the Flo- 
quet states and it illuminates the physical origin of the 
photon- assisted transport. Second method takes advan- 
tage of "Floquet Dyson's equation" to give an elegant so- 
lution which is more convenient for numerical solutioiP^. 
By taking the off-resonant limit of these solutions, the 
equivalence of transport properties under the application 



of light and those with effective photon-dresseed Hamil- 
tonian is established. 

Most of the analysis in this paper assumes the ab- 
sence of interactions among electrons as well as electron- 
phonon interactions. We argue in Section |IV| that, in 



the case of graphene and topological insulators under the 
off-resonant light, the results given in Section [IT] are ro- 
bust against these interaction effects at low temperatures. 
While the measurements of transport properties require 
the attachment of leads, the probe of the effective gap 
induced by light is plausible even in an isolated system. 
We propose in Section [V| such measurements through the 
adiabatic preparation of non-equilibrium systems com- 
bined with the transmission of probe laser with small 
frequencies. The essential ingredients in the arguments of 
Section [rVl and Section [Vl are the extensions of adiabatic 
theorem and Fermi golden rule to periodically driven sys- 
tems and Floquet states, dubbed as "Floquet adiabatic 
theorem" and "Floquet Fermi golden rule." We give the 
detailed proof of these important statements in the Ap- 
pendix. In Section [Vlj we conclude with possible exten- 
sions of this work. 



II. SUMMARY OF RESULTS 
A. Garphene effective Hamiltonian 

Here we consider graphene as an example of semi- 
metals and study the change in the transport properties 
under the application of light. We model graphene by a 
hexagonal tight-binding model with two 7r-bands, where 
we first neglect the electron-electron as well as electron- 
phonon interactions. In Section|lV|, we discuss the effects 
of these interactions and argue that they do not change 
the qualitative results of the analysis. We consider the 
application of circularly polarized light perpendicular to 
the plane of graphene. For concreteness, here we rep- 
resent the rotating electric field due to light as a time- 
dependent vector potential A(t) = 74(± sin(r^t), cos(r^t)) 
with E(t) = dA{t)/dt, where ft is the frequency of light. 
Plus sign is for right circulation of light and minus sign 
for left circulation. The light intensity is characterized 
by the dimensionless number A = eAa/h where e is the 
electron charge and a ~ 2. 46 A is the lattice constant of 
graphene. For intensity of lasers and pulses available in 
the frequency regime of our interests ~ lOOOTHz, A is 
typically less than 1. In this gauge, electrons accumulate 
phases as they hop in the lattice; 
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where Aij{t) = e/h{rj — r^) • A(t) with being the coor- 
dinates of the lattice site z, J is the hopping amplitude of 
electrons, and s =t, | are spins of electrons. For simplic- 
ity, we only consider the orbital effect of electromagnetic 
fields on electrons, and disregard the small Zeeman effect. 
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FIG. 2: a) The modification of tfie Hamiltonian due to tfie 
virtual pfioton process can be intuitively understood as the 
sum of two second order processes where electrons absorbs 
and then emits a photon, and electrons first emit and ab- 
sorb a photon, b) The illustration of the structure of ifeff 
in real space for graphene under the application of right cir- 
cularly polarized light in Eq. ([3|. The commutator [Hi, Hi] 
is the second neighbor hopping with phase cp = 7r/2. Thus, 
the tight-binding model under the application of light realizes 
Haldane model proposed in Ref. [7|. 



The inclusion of the Zeeman effect is straightforward. In 
this limit, spins trivially double the Hilbert space and 
thus we suppress the spin indices in the following. 

When the light frequency is off-resonant for any elec- 
tron transitions, light does not directly excite electrons 
and instead, effectively modifies the electron band struc- 
tures through virtual photon absorption processes. Such 
off-resonant condition is satisfied for the frequency Q ^ 
J in our model with 7r-bands. More general case of on- 
resonant light can be analyzed through the formalism 

The influence of 



developed in later section Section III 



such off-resonant light is captured in the static effective 
Hamiltonian i^eft^^ defined through the evolution opera- 
tor U of the system after one period T = 27r/r^ as 



^eff 



T 



\og{U) 



(2) 



where U = Texp ^— z H{t)dt^ and T is the time- 
ordering operator. Intuitively, i^eff describes the dynam- 
ics of the system on time scales much longer than T. In 
the limit of <C 1, i^eff is particularly simple near the 
Dirac points; 



^eff 



(3) 



n 



(for infinite system), 

where Hn is the discrete Fourier component of Hamilto- 
nian, i.e. Hn = ^ H{t)e'^^^^dt. In the second line, 
vg = 3 J/2 is the velocity of Dirac electrons, kx and ky 
are momenta measured from the Dirac points, cr^ and 
are Pauli matrices representing sublattice and valley de- 
grees of freedom, respectively. The modification of the 
Hamiltonian with respect to the static component Hq 



is the second term in Eq. (|3|. This term can be easily 
understood as the sum of two second order processes as 
illustrated in Fig. [2] a); one where electron absorbs a pho- 
ton and then emits a photon Hi iJ-i where co is 
the energy of the original electron; another where elec- 
trons first emits a photon and then absorbs a photon, 
which leads to H- 



-(uj-n) summing these two 

contributions, we obtain the correction due to the sec- 
ond order process, given in the second term of Eq. (|3|. 
In the second line, the plus sign is for right circulation 
of light polarization and the minus sign is for left cir- 
culation. For details of derivations, see Section IIIIl We 



note that the expression of the effective Hamiltonian in 
Eq. Q is only valid in the gauge in which light is repre- 
sented as time-dependent vector potential, and the effec- 
tive Hamiltonian has different forms for other gauge such 
as the one in which light represented as time-dependent 
electric fields. The effect of virtual photon absorptions at 
the degenerate Dirac points is to open a gap with mag- 



nitude A 



In Fig. Ji 



we illustrate the opening 



of the gap near one of the iJirac points upon the appli- 
cation of light for both infinite and finite systems. This 
Hamiltonian i^eff corresponds to a quantum Hall insula- 
tor, where each band is characterized by non-zero Chern 
numbei^^^. In Fig. [sjb), we have plotted the spectrum 
of i^eff where the system is infinite in x direction and 150 
sites in y direction with armchair edges. Here, we have 
chosen the intensity and frequency of light to be ^ = 0.3 
and Q = 7.5J. As a result of non-zero Chern number of 
the bands, i^eff shows the existence of gapless chiral edge 
states, colored as blue and green, corresponding to the 
edge state in the upper and lower edge, respectively. 

It is instructive to write the effective Hamiltonian in 
Eq. ([3| in real space. In the lowest order in A, Hi 
and H-i are the hopping between nearest neighbors with 
phase accumulations that depend on the direction of the 
hopping. Their commutators contain the second neigh- 
bor hopping with amplitudes q g^^/^ g^g illustrated 
in Fig. [2]b). Thus, the effective Hamiltonian is nothing 
but the Hamiltonian proposed by Haldane'' with sub- 
lattice potential M = and second hopping strength 



t2 



with the flux (f = 7r/2. 

above differ from that of Inoue and 
In their work, the 
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Our results 
Tanaka[18 in some important ways, 
effect of circularly polarized light on the Haldane model^ 
has been considered. They focused on the zero-photon 
sector of the Hamiltonian and concluded that the Chern 
number is zero whenever the second neighbor hopping 
t2 is zero or the staggered magnetic flux cp is zero, as is 
presented in Eq. (7) of their paper. Their work showed 
no transition from topologically trivial band insulators to 
topological non-trivial bands with Chern numbers. Here 
we considered a simple tight-binding model, correspond- 
ing to ^2 = and (p = of Haldane model. In con- 
trast with the result of Ref. [18], we show above (also 



see Section HI B 2 ) that the virtual photon absorption 



and emission process represented by the second term of 
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FIG. 3: a) The spectrum of graphene for a single spin near 
one of the Dirac points for an infinite system. Here we plot 
the energy spectrum of a static system corresponding to ^ = 
(upper figure) and the spectrum of iifeff for the system under 
the application of light with light intensity A = 0.3 (lower 
figure). A finite gap A opens at the Dirac point, b) The 
spectrum of ifeff for a single spin near one of the Dirac points 
for infinite system in x direction and 150 sites in y direction for 
A = 0.3. The frequency of light is chosen to be off-resonant 
with Q — 7.5 J. Here we chose the armchair edges, and in this 
case, the Dirac points are at = 0. ifeff shows the existence 
of gapless chiral edge states for each spin originating from 
non-zero Chern numbers of the bands, that are colored as 
blue and green, corresponding to the edge state in the upper 
and lower edge, respectively. The propagation of chiral edge 
states for a single spin is illustrated in the upper figure. 
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FIG. 4: Illustration of the response conduction current in 
Floquet Landauer formula given in Eq. The transmis- 

sion of the electrons can now happen with n photon absorp- 
tion/emissions. The total conductance is simply given by the 
sum of these contributions for each n. 



light as in Fig.[^. The non-equilibrium transport prop- 
erties of mesoscopic, periodically driven systems in this 
configuration have been studied previously, using Floquet 
theor 3^ combined with the Keldysh formalism^^. We ex- 
press the general result obtained in these works as the ex- 
tension of the multi-probe Biittiker-Landauer formulsP^ 
to periodically driven systems, a "Floquet Landauer for- 
mula" ; 



j: 



dc 



(4) 

(5) 



Eq. ([3| induces a gap at the Dirac point and leads to 
a non-zero Chern number. Such photon-dressing, which 
was neglected in the study of Ref . [18] , has dramatic ef- 
fects at the degenerate Dirac points and should be taken 
into account. 

There are a few different ways to probe the gap A in 
the effective Hamiltonian. For example, the gap opening 
might be confirmed through the observation of the trans- 
mission of low frequency probe lasers. In Section |V| we 
propose the possibility to probe this dynamically opened 
gap in an isolated system through adiabatic preparation 
of Floquet states. The main focus of this paper, how- 
ever, is the study of the manifestations of the gapped 
Hamiltonian i^eff in many-terminal transport measure- 
ments depicted in Fig. [l] which we now describe. 



B. Floquet-Landauer formalism and Hall current 

One of the central results of this paper is the demon- 
stration that measurements of DC current of the non- 
equilibrium system under the application of off-resonant 
light is determined by the static, photon-dressed Hamil- 
tonian i^eff- To this end, we consider the many-terminal 
measurements of DC current under the application of 



Here J^^ is the DC component of the current at lead a, 

= T /o Ja{t)dt. Here we assumed that the reservoirs 
are at zero temperature, and their chemical potentials 
{fia} are near the Dirac points, i.e. /j^a ^ 0. The DC 
current of a periodically driven system in Eq. Q con- 
sists of two physically distinct contributions; the pump 
current^^, J^^"^^, which can be present even when all the 
reservoirs have the same chemical potential /la and the 
response current, J^®®, which arises from the response 
of the driven system to the chemical potential differ- 
ences of the reservoirs. For the application of light we 
consider in this paper, the pump current J^^^^ is zero 
for inversion symmetric geometries, so we focus on the 
properties of the response current J^^ in the following. 
The transmission coefficients T^^ (n) in Eq. ([5| represent 
the transmission of electrons with energy ^ from 
lead P to a during which electrons absorb (emit) n pho- 
tons as illustrated in Fig. |4j Thus, the response con- 
duction in the presence of periodic drive can be seen as 
an extension of the static conduction, where the trans- 
mission can now happen with the absorptions/emissions 
of photons. We note that the expression of DC current 
in Eq. ([5| is valid for arbitrary strength and frequency 
of the drives. The transmission probabilities Ta(3{n) can 
be efficiently computed by dressing the propagators with 
photon absorptions/emissions-^^, and will be described 
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in Section |III[ We emphasize that the response current 
is generahy the sum of the contributions from n pho- 
ton absorption/emission processes (see Eq. ([5|) and thus 
its transport property cannot be described by any static 
Hamiltonians. The off-resonant case described below is 
an except ionahy simple case in this respect. 

We employ this Floquet-Landauer formalism to study 
the off-resonant, large frequency regime J <^ Q with 
weak intensity of light, i.e. J\? <C 1. In this regime, 
absorptions or emissions of photons are suppressed by 
A^^ and the transmission coefficients Tai3{n) with n ^ 
is small and of the order of 0{A'^). On the other hand, 
the zero-photon absorption/emission transmission coeffi- 
cient Taf3{0) is modified due to virtual photon processes. 
Such modifications are included in i^eff and the transmis- 
sion probability is given by Tq,^(0) = T^^ + 0(^^), where 
T^^ is the transmission probability of the static system 
described by i^eff- These results will be rigorously estab- 
lished in Section Hill 

This correspondence demonstrates, under our assump- 
tions, that graphene under the application of off-resonant 

light behaves as an insulator with gap A = ^^^^ with 
Hall conductance quantized at 2e^ /h with possible cor- 
rections up to the order of 0{A^). Here the factor of 
2 comes from spin degrees of freedom. While we estab- 
lished the results in the perturbation theory on A^ it 
is possible to analytically confirm the insulating behav- 
ior for all orders in A for weak contact couplings with 
leads (see Section HI). We emphasize that although the 
effective Hamiltonian is perturbative in A^ the Hall con- 
ductance at zero temperature is non-pert urbative: an in- 
finitesimal gap A is sufficient to yield a topological band 
with non-zero Chern number. 

A distinct feature of this light-induced Hall effect above 
is that the Hall conductance switches its sign under the 
change of circulations of light polarization. This can be 
easily checked for the geometry of the system which is 
symmetric under x ^ — x, under which the circulation of 
light reverses. Such reversal of Hall current can be used in 
the experiments to distinguish this light-polarization de- 
pendent current from light-polarization independent cur- 
rent, which could originate from mechanisms we did not 
consider in this paper. 

We briefly describe the requirements to observe the 
proposed phenomena with off-resonant light in graphene. 
The band width of graphene in the tt orbital is given 
by 6 J where J ^ 2.4eV, placing the required frequency 
of off-resonant light to be soft x-ray regime with Q = 
3500THz. For this frequency of light, the gap of the sys- 
tem A can reach A ^ 300K for the strong light intensity 
I ^ 3x lO^^W/cm^^^ which gives A ^ 0.09, where we ex- 
pect the Hall conductance to be quantized with possible 
correction of 1% of 2e^ /h. In reality, even such high fre- 
quency of light is expected to be absorbed in graphene. 
Such direct electron excitations lead to reconfiguration 
of electron occupation numbers, which modifies the Hall 
conductance from its quantized values. 



C. Three dimensional topological insulators 

The analysis of graphene above can be directly ex- 
tended to three dimensional topological insulators such 
as Bi2Se3. The low energy description of electrons on 
surfaces of Bi2Se3 is given by two dimensional Dirac 
fermion^^, and is described by the Hamiltonian H^^"^^ = 
VTiikxC^y — kyax) where vti is the velocity of the Dirac 
fermion, and are Pauli matrices corresponding to two 
bands near the Dirac point. As before, we assume the 
application of weak, off-resonant, circularly polarized 
light. The orbital effect of the light is taken into ac- 
count through the replacement k ^ k — A(t). At the 
Dirac cone, the virtual photon process again opens a gap 
and the effective Hamiltonian is (see Eq. (Is])) 



surf 



eff 



A^v^ 



TT 



(6) 



where +(— ) corresponds to the gap due to right (left) 
circularly polarized light. The consequences of the gap 
coming from the third term in Eq. ([6| are extensively in- 
vestigated in Ref [8 21 . Just as in the case of graphene, 
the induced insulator is topologically non-trivial, and ex- 
pected to result in anomalous quantum Hall effect with 

2 

Hall conductance with possible corrections up to the 
order of 0{A^). Here we propose to probe the unique 
magneto-electric response of the gapped topological in- 
sulator through pump-probe type measurements, where 
circularly polarized light is used to open a gap at the 
Dirac point and linearly polarized light with small fre- 
quency within the gap is used to probe the Faraday /Kerr 
rotations^^, as illustrated in Fig.jsja) (also see Section[v|). 
Unlike other schemes proposed previously with ferromag- 
netic layers, here the Faraday/Kerr rotations can only re- 
sult from the topological insulators and they give unam- 
biguous signature of magneto-electric effects. In a simi- 
lar fashion, the existence of magnetic monopoles can be 
probed by placing an electric charge near the surface of 
the topological insulator in the presence of circularly po- 
larized light (see Fig. |5]b)). 



D. Discussion 

In the analysis of graphene and topological insulators 
above, we assumed the off-resonance of light for entire 
bands, but the gap in the effective Hamiltonian opens 
whenever the light is off-resonant near the Dirac points, 
which requires much less stringent condition on the light 
frequency. However when inter-band electronic transi- 
tions occur due to photon absorptions, subsequent relax- 
ation processes are expected to change the electron occu- 
pation numbers in the steady state, and modify Hall con- 
ductance away from quantized values. Thus in the case of 
on-resonant light, the system is expected to display non- 
quantized Hall effects without magnetic fields. Moreover, 
as we show in Section llllj the application of on-resonant 
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FIG. 5: a) Measurements of Faraday/Kerr rotations in three 
dimensional topological insulator. Circularly polarized light 
with a large frequency is used to open a gap at the Dirac 
point and linearly polarized light with small frequency within 
the gap is used to measure the Faraday/Kerr rotations. The 
polarization angle of the light is denoted by blue arrows in 
the figure, b) The induction of magnetic monopole inside a 
three dimensional topological insulator through the applica- 
tion of light. Electron charge is placed near the surface, and 
circularly polarized light is applied on the surface to break the 
time-reversal symmetry and open the effective gap. The mag- 
netic monopole is induced as a mirror image of the electron 
charge. 



light leads to the photo-assisted conductance and the re- 
sulting non-equilibrium transport property can no longer 
simply be described by the static effective Hamiltonian. 
The transport under the on-resonant light contains rich 
physics in itself, and it will be studied in the future 
works. On the other hand, it is possible to achieve the off- 
resonance with small frequency of light in, for example, 
the gapped systems such as Boron-Nitride by applying 
the sub-gap frequency of light. Many of these possible 
extensions can be studied through the formalism devel- 
oped in the following sections. The understanding and 
formalism obtained in this paper is likely to guide future 
searches for the optimal systems to study photo-driven 
quantum Hall effects without magnetic field. 

III. NON-EQUILIBRIUM TRANSPORT: 
FORMALISM 



where ja are the sites in graphene that are connected 
with leads a, and Ta{uj) represents the coupling strength 
with leads Ta{cjo) = tl^pa{oo) where ta is the hopping 
strength from graphene to lead a and pa{^) is the density 
of states in the lead a at energy uj. Also, fa{^) is the 
Fermi function at lead a, fa{^) = g/3c.(u;-/xc.)+i where 
Pa = ^/ksTa is the inverse temperature and pa is the 
chemical potential of the reservoir connected to lead a. 
If we take the zero temperature limit and assume that 
differences of chemical potentials at each lead are small, 
the expression in Eq. ^ is reduced to the simpler Floquet 
Landauer formula given in Eq. ([5|. 

The calculation of conductance given by J^^ reduces 
to the calculation of Floquet Green's function Q{uj,n). 
Here, Qi^i^ (a;, n) is nothing but a Fourier transform of the 
retarded Green's function Gf'i>{t^t'). Starting from the 
usual definition of the retarded Green's function, 

Gf,(t,i') = -ie{t-t'){{ci{t)c\,{t')) + {c\,{t')ci{t))) (9) 
we take the Fourier transform to obtain 

/oo 
dt'G/fp(t,t')e'^"+'''^^<*-*'^- (10) 
-OO 

Because we are driving the system at the given frequency 
1], this Green's function, as a function of t, should contain 
only the discrete Fourier components. Therefore, we can 
expand as 



CO 

G5,(t,c^)= Gui'Me-'^''' (11) 

/c= — OO 



A. Floquet Landauer formula 

In this paper, we study the transport properties of sys- 
tems under the application light in the Landauer-type 
configuration, where the systems are attached to the 
leads as in Fig. [l] Previous works^^ obtained the DC 
current in periodically driven systems in terms of the 
Floquet Green's functions Q{uj,n) which represent the 
Fourier transform of retarded Green's function, and is 
a propagator with frequency uj which absorbs (emits) n 
photons. Here we express the general results in physically 
transparent fashion; 

jdc rpump , rres 

= EE / ^T^An,u:)iUH-Uuj)) (7) 



The equation of motion followed by Q{uj^n) can be 
obtained by writing out the equation of motion for 
Gf'i,{t^t') and taking its Fourier transform. Th result- 
ing equation can be written in the most compact form in 
the matrix equation whose elements correspond to differ- 
ent (discrete) frequency components, nVt. Explicitly, the 
equation is given by 



(cj + l^-H-ir/2)g = I (12) 



where 
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(13) 



f(w + n) 



f(w - n) 



/ . 



5(0, w) g{i,Lo) 
g{-i,uj) g{o,uj) g{i,uj) 
g{-i,uj) a(o,w) 



■ / 



H is the matrix of Hamiltonian, whose elements are 
the Fourier components of Hamiltonian given by Hn = 
^ H{t)e^^^^dt. Q, is the diagonal matrix whose el- 
ement is just the discrete frequency of the component, 
r is a diagonal matrix that represents coupling of the 
systems with leads. Its element r(co') has non-zero value 
only at site that couples with leads, and with value 
= ^cx{^) = taPa{^)' Finally, g is the matrix 
composed of the Floquet Green's functions. 



The equation in Eq. (12) can be thought of as the ex- 



tension of the equation of motion for free electrons cou- 
pled with leads to periodically driven systems. The rest 
of this section is devoted to the solution of this equation, 
and to the explanations of its physical significance for the 
response current given by Eq. 



B. Floquet states and Floquet Green's functions 

1. General relation 

The "stationary states" of the Schrodinger equation 
for periodically driven systems are the states which re- 
turn to themselves after one period of time, T = 27r/l], 
with possible phase accumulations. These so-called Flo- 
quet states are the eigenstates of the evolution operator 
over one period, and thus also eigenstates of effective 
Hamiltonian i^eff defined in Eq. ([2|. Green's functions 
Q{uj^n) that describes the propagation of particles with 
possible absorptions/emissions of photons have natural 
expressions in terms of these Floquet states. 

The time evolution of the Floquet states \if>a{f)) can 
be expanded in the discrete Fourier component of the 
driving frequency ^7, and can be expressed as 



\^a(t)) = e 



-iEat 



-iQnt 



(14) 



In the following, we give two different solutions of 
Eq. (12). In Section IIIB, we solve the equation by 



expressing the Green's functions in terms of Floquet 
states, the "stationary states" of periodically driven sys- 
tems after one period of time. This solution illustrates 
the physical origin of the transport given in Eq. (0. 
In Section IIIB 2[ we derive the equivalence of non- 
equilibrium transport and the transport given by the ef- 
fective photon-dressed Hamiltonian i^eff claimed in Sec- 
tion [H] by taking the off-resonant and weak intensity 
limit. In Section III C[ we give another solution of 



Eq. ( 12 ), which is valid for a certain class of periodic drive 



including the application of circularly polarized light. 
This solution is derived by writing "Floquet Dyson's 
equation," and has the advantage of being numerically 
efficient. Again, by taking the off-resonant, weak inten- 
sity limit of this solution, we arrive the result reported 
in Section HTl 



where Ea is the quasi-energy of the Floquet state \(pa) for 
i^eff and \(p'2) is nth Fourier component of the Floquet 
state. The eigenstates of effective Hamiltonian i^eff are 
given as \(pa) = X^n l^a)- can see from the ex- 

pression above, the quasi-energy Ea is only well-defined 
up to the driving frequency 1^, i.e we can equally de- 
fine Ea-\-mQ as the quasi-energy of \(Pa{t)) by redefining 
l^a) ^ \^a~^)' Physically, this means the quasi-energy 
is only conserved up to the driving frequency Q because 
the system can absorb or emit photon energies. Also, this 
fact can be seen as a natural consequence of the break- 
ing of continuous time-translation invariance through ex- 
ternal drivings where the system only possesses discrete 
time-translation invariance under t ^ t -\- T. In the fol- 
lowing, we assume that —Q/2 < Ea < without loss 
of generality. We take the normalization of Floquet states 
such that Xln(^al^6) — ^ab- The Schrodinger equation 
for the Fourier components of the Floquet states is time- 
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independent, 

{E, + n^l)\^l)=Y,H^.mm (15) 

771 

where is the discrete Fourier component of Hamil- 
tonian, i.e. = ^ H{t)e^^^^dt. This equation en- 
capsulates the evolution of states that allows the absorp- 
tions/emissions of photons; application of Hamiltonian 
leads to the absorption of m photons and the state 
\ip'^) is the component of the state with n photons. Here 
we considered the evolution of the systems in the absence 
of coupling with leads, but in the presence of the coupling 
with leads, zero frequency component of the Hamiltonian 
Hq contains the imaginary "leaking" term iV /2. 

This Schrodinger equation takes, in the matrix form. 



(16) 



where 



V : J 



(17) 



The Hamiltonian matrix H and driving frequency matrix 
are given in Eq. (13). 



Thus, Floquet state \ipa) is nothing but the eigenstates 
of the composite Hamiltonian (—11 + H). Notice that 
"shifted" state 



\ 



V 



(18) 



is also an eigenstate with eigenvalue Ea + nVt. 

Now we can relate the Floquet eigenstates given by 
Eq. (16) and Floquet Green's function given by Eq. (12). 
The formal solution of Eq. (12) is obtained in terms of 
the eigenstates \Lpa) of the matrix 11 — H — ir/2 as 



= E 



l</'a)(<,5a| 



(19) 



Here, {(pa\ is the state that is determined from {(pa\^h) = 
5 ah' (Note that in the presence of ir/2, {(f>a\ is not just 
a complex conjugate of As is clear from Eq. (16), 

the eigenstates of \ifa) of the matrix 11 — H — ir/2is 
nothing but the Floquet states in the presence of the 
coupling with leads, represented by —iT/2. Notice that 
the eigenstates \ipa) in Eq. (19) include all the shifted 



understanding, we obtain the expression of the Floquet 
Green's function Q{n^uj) as 



Ea — mVt 



(20) 



states |(/9^^^^*(n)) in Eq. (18) for all integers n. With this 



We can see that Floquet Green's function is an intuitive 
extension of the Green's function for free electrons that 
allows absorptions and emissions of photons. As is ex- 
pected, this Green's function transfers the state in —m 
photon sector {(Pa^\ to n — m photon sector \(Pa~^) by 
absorbing n photons. 

In the presence of on-resonant light, Floquet states 
generally contain non-zero amplitudes in \(f^) for more 
than one value of n, and therefore, the contributions to 
the response current in Eq. ([5| from a few photon ab- 
sorptions/emissions are non-zero. Thus effective static 
Hamiltonian or band structures, which are only the de- 
scription of average of n photon states \^a)^ does not 
appropriately capture transport properties under the on- 
resonant light. In this case, it is necessary to compute 
the full Floquet Green's function given by Eq. (20) and 
calculate the response current in Eq. 

On the other hand, in the case of off-resonant, weak 
intensity of light, transport properties of non-equilibrium 
systems can be described by an effective photon-dressed 
Hamiltonian. In the next section, we provide the proof in 
the case of semi-metals such as graphene and topological 
insulators. 



2. Effective Hamiltonian description 

As summarized in Section [TTj a rich physics appears 
when circularly polarized light is applied to graphene 
and topological insulators. The description of the non- 
equilibrium transport takes a particularly simple form 
for the off-resonant light in the limit of small light in- 
tensity ^ <C 1. Here we apply the general formalism 
developed in the previous section to these systems and 
study the transport property by obtaining the Floquet 
Green's function in Eq. ( [2Q| ). 

From the explicit form of the Hamiltonian in Eq. ([T]), 
it is clear that Hn ~ 0(^1^'). This simply means 
the absorptions of photons are suppressed by the factor 
0{A^^^). Thus, the n photon sectors of the Floquet states 
are expected to scale as \(f2) ~ 0(^1^') with zeroth 
order solution being the static part of the Hamiltonian 
Hq -\- ir/2. Here the term ir/2 represents the coupling 
with leads and we assumed the same strength of the cou- 
pling r = Fc^ at each lead a and further assumed that 
it is independent of frequency. This latter assumption 
is not important in the off-resonant case because cur- 
rent is essentially conducted only at chemical potential 
of leads, as we will confirm later. Starting from the equa- 
tion Eq. (15), we apply a degenerate perturbation theory 
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in the lowest non-trivial order in A to obtain 



1 



nf2 



(21) 
(22) 



In the derivation, we assumed <C 1^, so the expression 
above is only valid near the Dirac points. Note that since 
the Hamiltonian Hq is degenerate at the Dirac points, the 
mixings of the states due to the perturbations of A are 
not small. This result indeed shows that \(fa)^ 0(^1"^'), 
and therefore, the Floquet states \(pa) can be approx- 
imated by the zeroth level of the Floquet states \(Pa)^ 
which is given by the eigenstates of the effective Hamil- 
tonian i^eff = + — — plus the coupling with leads 

ir/2. 

Using the solution of Floquet states above, we can ob- 
tain the response current in the lowest order in A. The 
scaling |(/?^) ~ 0{A^) in Eq. (22) directly implies that 
|^(n,cj)p ~ 0{A'^^). Moreover, Green's function with 
no photon absorptions or emissions can be approximated 



Y^\M§1^0{A') 



Ea 

h 0{A^) 



where \ip^^) is the eigenstates of i^eff + if /2 , and there- 
fore Q^^{(jj) is the free electron Green's function for the 
static system with Hamiltonian i^eff coupled with leads. 
Thus, these arguments combined with the expressions 
of currents in Eq. ^ and Eq. ([8| prove that the many- 
terminal measurements of non-equilibrium systems under 
the off-resonant light give the same result as if the system 
is given by the static effective Hamiltonian i^eff- We em- 
phasize that this result immediately implies the following 
two facts; 1. the non-equilibrium system displays insu- 
lating behaviors in the longitudinal conductance 2. Hall 
conductance is nearly quantized with possible correction 
of 0{A?)^ as presented in Section [h| 



3. Insulating behavior for gapped effective Hamiltonian Hejj 

In the analysis of the previous section, we established 
the insulating behaviors of graphene under the applica- 
tion of off-resonant light through the perturbation the- 
ory in A. Such analysis only shows that the longitudinal 
conductivity is small and of the order of but does 
not show, in a strict sense, that the conductivity goes 
to zero at zero temperature. Using the formalism devel- 
oped in previous sections, it is possible to show that the 
non-equilibrium system is an insulator for all orders in 
A as long as the chemical potential of leads lies below 
the effective gap of i^eff- The argument does not rely 



> 



C/5 



> 
I 



C/5 



FIG. 6: Illustration of the configuration considered in Sec- 
tion jlll iTSl to demonstrate the insulating behavior for gapped 
effective Hamiltonian i/eff- Nl number of leads are attached 
to the infinite plane of graphene as "left" leads and Nr num- 
ber of leads are attached as "right" leads. The chemical po- 
tentials of the reservoirs connected to left leads are assumed 
to be at V/2 and those of the reservoirs connected to the right 
leads are at -V/2 with \V\ <^ J 



on the off-resonant condition and in principle applicable 
whenever i^eff has a gap. 

Here we consider an infinite plane of graphene, and we 
attach Nl number of leads as "left" leads and Nn number 
of leads as "right" leads, where these leads are separated 
by a large distance, see Fig. |6] Here we assume that the 
leads are coupled with the system with equal strength, 
given by r(co'). For clarity, we consider the situation in 
which the chemical potentials of the reservoirs connected 
to left leads are at V/2 and those of the reservoirs con- 
nected to the right leads are at —V/2 with |V| ^ J. 

In the limit of small coupling strength, the Green's 
function in Eq. (20) can be obtained through the pertur- 



bation theory on F, and are given by 



(23) 



where 7a(w) = E„«|f (c^+nl^)/2|0, and |0 and^^ 
are the Floquet states and (quasi-) energy of the system 
in the absence of the coupling with leads. In the limit of 
small 7a (^), the square of the Green's function Q{n^uj) 
can be approximated by a delta function, so that the 
transmission probability also becomes a delta function in 
frequency; 



J{3LM-n\'\{^a"'\jR.)\' 



7tS{uj — Ea — mQ) 



27a (C^) 

where jLi{Ri) are the sites of left (right) leads. Now 
note that the (quasi-) energies Ea are the eigenvalues 
of i^eff in Eq. (2) in the main text. Therefore, if all 
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the chemical potentials lies within the gap of i^eff^ ^•6- 
\V\ < A the delta function gives zero everywhere for 
— V/2 < uj < V/2. Note that we have taken the quasi- 
energies Ea to lie between —Q/2 < and thus, 

by assumption, A < Q. mQ term in the delta func- 
tion of Eq. (24) accounts for the possible transmission 



< 



+ 



of electrons at high/low energies through photon absorp- 
tion/emission processes. As long as we are interested in 
the transmission of electrons near the chemical potential 
which lies within the effective gap, mQ term plays no role 
in the conduction. Thus, from the expression of the DC 
current in Eq. ([7|), it is clear that the current has to be 
zero when |V| < A. 

The argument above is general and did not require the 
condition of off-resonance. In the case of on-resonant 
light, the effective band structures are given by mixing 
the static eigenstates whose energies differ by This 
"folding" of the band structures generically leads to a 
large number of states appearing in the effective Hamil- 
tonian near the chemical potential, and subsequently the 
effective gap in i^eff becomes proportional to 0{A^) with 
n approximately determined by the ratio of the static 
band width and the driving frequency Q. While the insu- 
lating behavior should be observable in the small window 
of the gap, the gap could be small in this case. 



C. Floquet Dyson's equation 

In this section, we present yet another way to obtain 
the Floquet Green's function Q{n^uo) which gives an ef- 
ficient way to numerically evaluate the Floquet Green's 
function for a certain class of periodical drives. 

In this section, we consider the Hamiltonian that de- 
pends only on the first harmonics of the driving frequency 
r^, namely, the Hamiltonian takes the form 



H{t) = Ho^ Vie- 



-iQt 



VLie 



(24) 



For example, for the application of the circularly po- 
larized light to two dimensional lattice systems, Vi = 

"^Zji^j + V-i = Vi in the gauge in which 

the light is represented as a circulating potential. How- 
ever, in this gauge, V±i diverges as Xj^yj oo, care must 
be taken to study with this gauge. Conceptually useful 
gauge is the gauge in which the effect of light is repre- 
sented as a phase accumulation as in Section |n[ For weak 
amplitude of light, we can approximate the Hamiltonian 
in this gauge in the form of Eq. ( 24 ) with Vi = Hi and 
Vi = H_i. As before, we are interested in the terminal 
measurements of conductance, and thus assume that the 
static part of the Hamiltonian Hq contains the "leaking" 
of particles into leads given by iV{uo)/2. 

In order evaluate Floquet Green's functions, we first 
rewrite the equation Eq. (12) in the suggestive form of 



V-i 



< = + < = 



FIG. 7: The Floquet Dyson equation. The propagator goes 
from right to left. Double line represents the full propaga- 
tor G{n,(jj) and single line is a bare propagator G{n,(jj) = 
uj-\-nn-Ho which does not include the effect of photon absorp- 
tions. 



"Floquet Dyson's equation" (see Fig. [t]); 

G°(n, n) (V-iGin + 1, c^) + ViGin - 1, cj)) 

(25) 

Here G^(n,l]) = represents the bare propa- 

gator of a particle with n photons. This equation has 
the intuitive understanding of the full propagator ^(n, uj) 
that represents the n photon absorption process as being 
composed of the full propagation of ^(n ± l,cj) followed 
by the absorption or emission of a photon, followed by 
the propagation of the bare particle. 

A particularly elegant solution for G{n^uj) is provided 
by continued fraction method^^. The building block of 
the solution is the dressed propagator 



J^_(n, uj) 



1 



{G^)-\n,uj)-V-i 



1 Y 

(GO)-i(n+l,a;)-y_i:^yi 

for n > 

1 



for n < 

The propagator JFj^{n^uj) is dressed only from the higher 
photon number states, and the propagator J^-{n^uj) is 
dressed by the lower photon number states. The full 
propagator is then given as 

G{0,uj) = {uj-Ho-Vef,)-^ (26) 

Kff = ViA(-l,c^)F_i + F_iJ-+(l,c^)Fi 
G{n,uj) = :F+{n,uj)Vi'":F+{l,uj)ViG{0,uj) 

for n > 

= T-{n,uj)V-i ■ ■ ■ T-{-l,uj)V-iG{0,uj) 

for n < 

(27) 

This solution is valid for any driving frequency. Remark- 
ably, we see that the zero-photon absorption propagator 
G{0^uj) is simply given by the propagator in an effective 
Hamiltonian i^eff = + ^ff- 
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For the gauge in which hght is represented as time- 
dependent vector potential, and weak intensity of hght 
^ <C 1, we can approximate J^+(l,cj) = l,a;) = 

^ in the hmit of high frequency. Thus, we reproduce 
the result we obtained in Section [111 B 21 of the effective 

\H H 1 

Hamiltonian i^eff = + ~q — — in this limit. 



IV. EFFECT OF ELECTRON-ELECTRON AND 
ELECTRON-PHONON INTERACTIONS 

The non-equilibrium transport properties described 
in Section |TT] are robust against interactions such as 
electron-electron interactions, interactions between elec- 
trons and disorder, and electron-phonon interactions. 
The electron-electron interactions only renormalize the 
velocity of Dirac electrons, vq and vjt, and do not change 
the Dirac nature of the electrons near the Fermi surface^^. 
The quantum Hall insulators are insensitive to disorders 
due to the topological origin of the phase, as long as the 
disorder strength is small compared to the gap size, A^^. 

The robustness of the phenomena against phonon scat- 
terings originates from the conservation of energy in i^eff 
up to the light frequency Vt. When the chemical poten- 
tials of leads lie in the gap of i^eff, the non-equilibrium 
current in many-terminal measurements is conducted 
through electrons in the lower band of i^eff- Such cur- 
rent can degrade due to electron-phonon interactions if 
electrons in the lower band can be excited to the higher 
band. However, such excitations in the bands of effec- 
tive Hamiltonians require a physical energy greater than 
the gap A as is rigorously established in "Floquet Fermi 
golden rule" in Appendix A. It is in principle possible 
to absorb energies from photons, but because the fre- 
quency of photons Vt is assumed to be much larger than 
band width, the absorption of such large energy requires 
the excitations of electrons together with many phonons, 
and therefore such a process is suppressed. Thus, the 
transition of an electron from the lower band of the ef- 
fective Hamiltonian to the higher band is possible only 
through the absorption of phonon energies. Therefore, at 
low temperatures, the property of an "insulating" state 
of the effective Hamiltonian is protected against electron- 
phonon interactions by the gap. 



V. PROBE OF THE INDUCED EFFECTIVE 
GAP IN AN ISOLATED SYSTEM 

In this section, we propose a different way to probe 
the effective gap induced by off-resonant light through 
the pump and probe measurements in an isolated sys- 
tem. The essential idea is simple. Given a system under 
the application of light (called pump laser), suppose that 
the effective Hamiltonian i^eff defined by Eq. Q has a 
gap A. We prepare the state, in isolation from thermal 
reservoirs, such that only the lower band of i^eff is occu- 
pied through a sort of "adiabatic preparation." Here we 



start from zero-temperature static system, and increase 
the intensity of light gradually to increase the size of the 
gap, A. As we argue below and Appendix A, adiabatic 
theorem in Floquet picture guarantees that the final state 
has the electron occupations such that the lower band of 
photon-dressed Hamiltonian i^eff is occupied. Now for 
this occupation of electrons with a gap to a higher band, 
it is intuitively clear that the system becomes transpar- 
ent to the probe light with frequency smaller than A. 
In the case of graphene and topological insulators under 
the application of light, we expect that the transmitted 
probe light results in the Faraday rotations. 

The pump and probe measurements described above 
are well-understood if the modification of the system 
from the original Hamiltonian to final Hamiltonian i^eff 
is done through a static fields. In this case, the adiabatic 
preparation is guaranteed by adiabatic theorem, and the 
transmission of probe light can be confirmed by looking 
at the Fermi golden rule which shows that the photons 
cannot be absorbed by electrons due to the conservation 
of energy. 

In the case of periodically varying fields, analogous 
statements hold. "Floquet adiabatic theorem" shows 
that, under an adiabatic evolution of the periodically 
varying fields, each Floquet state follows the instanta- 
neous Floquet state given by the instantaneous Hamil- 
tonian. Similarly, "Floquet Fermi golden rule" gives the 
rate in which the transition from one Floquet state to 
another happens under small perturbations. This result 
shows that the quasi-energies of Floquet states are con- 
served up to integer multiples of driving frequency Vt. 
Thus as we have claimed above, the electrons in the lower 
band of i^eff cannot be excited to higher bands unless the 
photon energy is larger than the band gap A. We give 
the detailed proof of these theorems in the Appendix A. 



VI. CONCLUSION 

In this paper, we studied the transport properties of 
non-equilibrium systems under the application of light 
in many-terminal measurements. Starting from Floquet- 
Landauer formula, we gave two different solutions of 
Floquet Green's functions that illustrate the physical 
origin of transport in this situation. We found that 
for generic driving frequencies, the transport involves 
photon-assisted conductance and cannot be described by 
any static, effective Hamiltonians. 

In the case of graphene and topological insulators un- 
der the off-resonant light, the non-equilibrium transport 
does not involve photon absorptions/emissions. Rather, 
the electron band structures are modified through the 
virtual photon absorption/emission processes. We es- 
tablished, through the solution of Floquet Green's func- 
tion, that such modifications are captured by the static 
photon-dressed Hamiltonian, and that the transport in 
this system becomes equivalent to that described by the 
photon-dressed Hamiltonian. Remarkably, the effective 
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Hamiltonian obtained in this way takes the form of Hal- 
dane modeF with second neighbor hopping with phase 
accumulations for graphene under the apphcation of cir- 
cularly polarized light. 

One important aspect of our proposal is the opening 
of the gap in the photon-dressed Hamiltonian when the 
original static Hamiltonian is semi-metal and gapless. We 
gave two physical manifestations of such a gap. One is 
the insulating be havior of the driven system attached to 
the leads (Section HI B 3 ). The attachment of leads is cru- 
cial to determine the electron occupation numbers. An- 
other is the transmission of low frequency light in an iso- 
lated sy stem (Sect ion |V|) after the adiabatic preparations 
of states. We argued the possibility of such pump-probe 
measurements by establishing two important extensions 
of well-known theorems, "Floquet adiabatic theorem" an 
"Floquet Fermi golden rule" (Appendix A). 

The formalism and intuitive understanding developed 
in this paper can be used to study the transport proper- 
ties of a variety of systems under the application of light. 
It is of interests to analyze, for example, the transport 
properties of light-induced topological systems proposed 
in Ref. [10]. In addition, our analysis shows that trans- 
port under the application of light contains richer physics 
than static transport. In particular, photon- assisted con- 
ductance in which electrons absorbs/emits photons dur- 
ing the propagations is the unique feature of driven sys- 
tems, and it is interesting to analyze how such physical 
process results in energy conductions. While we focused 
on the response current J^es in this paper, yet another 
aspect of driven systems is the presence of pump current 
-^pump appearing in Eq. Q . It is of interests to find ma- 
terials that can pump currents by simply shining light on 
their surface. 
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an analogous theorem as Fermi golden rule in the period- 
ically driven systems. When the chemical potentials of 
leads lie in the gap of i^eff, the non-equilibrium current in 
many-terminal measurements is conducted through elec- 
trons in the lower band of i^eff- Such current can de- 
grade due to electron-phonon interactions if electrons in 
the lower band can be excited to the higher band. By 
deriving "Floquet Fermi golden rule," we demonstrate 
that such excitations in the bands of effective Hamilto- 
nians still require a physical energy greater than the gap 
A. It is in principle possible to absorb energies from 
photons, but because the frequency of photons Q is as- 
sumed to be much larger than band width, the absorp- 
tion of such large energy requires the excitations of elec- 
trons and many phonons, and therefore such a process 
is suppressed. Thus, the transition of an electron from 
the lower band of the effective Hamiltonian to the higher 
band is possible only through the absorption of phonon 
energies. Therefore, at low temperatures, the property of 
an "insulating" state of the effective Hamiltonian is pro- 
tected against electron-phonon interactions by the gap. 

The proof of the second statement requires two steps. 
If we assume that the closed system with i^eff can be 
prepared in a state such that only the lower band of i^eff 
is occupied, then we can argue from the "Floquet Fermi 
golden rule" that the low frequency laser with uj < A 
cannot be absorbed by the electrons. Therefore, such a 
system is transparent to the light. In order to prepare 
such a "filled" state of the effective Hamiltonian i^eff, 
we consider an adiabatic preparation. Starting from the 
half-filled state of original systems whose chemical po- 
tential lies at the Dirac points, we adiabatically increase 
the strength of light. We argue, by deriving "Floquet 
adiabatic theorem"'^, that such procedure prepares the 
filled state of i^eff except possibly exactly at the Dirac 
points. 

These two statements rely on two general theorems 
about periodically driven systems, dubbed as "Floquet 
Fermi golden rule" and "Floquet adiabatic theorem." In 
the following, we derive these results, using the elegant 
approach from "two-time" formalism^^. We emphasize 
that these results are general and have wide applications 
outside of what we discussed in this paper. 



Appendix A: Floquet Fermi golden rule and Floquet 
adiabatic theorem 

In this section, we establish the following two state- 
ments about i^eff studied in the main text; 1) The result 
of many-terminal measurements of the systems under the 
application of light obtained in the main text is robust 
against electron-phonon interactions, as long as the en- 
ergy of phonons dictated by the temperature of the sys- 
tems is smaller than the induced gap A; 2) The photo- 
induced gap A can be probed, in a closed system, by the 
transmission of a laser with frequency cu < A. 

We give arguments for the first statement by deriving 



1. Two-time Schrodinger equation 

In order to study the dynamics of periodically driven 
systems, it is convenient to separate two time scales, a 
fast time scale associated with the driving frequency Q 
and a slow time scale associated with other dynamics 
such as those of phonons. We let t denote the former 
time scale and r the latter, and obtain the Schrodinger 
equation of the slower dynamics in terms of r through 
the replacement id/dt id/dt -\- id/dr. Then the time 
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evolution of states for slow time scale can be written as 

i^mr)) = (n + V{T))\4>{T)) (Al) 
d 



n = H{t)-i 



dt 



(A2) 



where H{t) corresponds to the Hamiltonian with peri- 
odic drives with frequency Vt and V{t) represents the 
perturbation of the system with slow frequencies com- 
pared to Vt. In the absence of the perturbation F(r), the 
eigenstates of the Schrodinger equation above is given by 
Floquet states such that 



(A3) 



where = e^^"^|(p«(t)) = En^~'^"'\^a) is a state 

with a periodic structure \^a{t)) = |^Q;(t + T)), and 
Ea is the quasi-energy of the Floquet state, i.e. the 
eigenenergy of i^eff- Here \(pa{t)) represents a Floquet 
state which satisfies the equation 1-L\(pa{t)) = 0. Note 
that Ea is only defined up to so that physically 
the same Floquet state \(pa{t)) in Eq. ([T4| can be as- 
sociated with the eigenvalue E^ + mQ and the state 
\^a{t)) = Ene"'^''%S^'^). Here we take the con- 
vention that \^'^{t)) with m = is associated with the 
quasi-energy E^ such that —Q/2< E^ < 0^/2. The or- 
thogonality of the eigenstates \^oc) can be recovered by 
defining the inner product of Floquet states as the aver- 
age of the usual inner product over one period of time, 



iiXaM) = ^ {Xa{t)\x^{t))dt. 



(A4) 



Then we have {{^a{t)\^^{t))) = 5^,^. 

These extensions of the inner products and eigenvalue 
problem in periodically driven systems can be consid- 
ered as the extension of Hilbert space to include the fast 
time variable t as another "spatial" variable. The inner 
product Eq. (A4) in this Hilbert space integrates over 



t, and the time variable is now represented by the slow 
time variable r. We point out that the inner product 



Eq. (A4) makes sense only when any dynamics associ- 



ated with r occurs in a slow time scale than the period 
of driving T. In principle, operators and states which 
depend on r change during the integration time T of fast 
time variable t due to the dependence on r. Since we 
are treating t an d r as independent variables, the inner 
product Eq. (A4) ignores such r dependence. As long as 



such changes are small, the inner product Eq. (A4) gives 
a good approximation. 

The crucial observation is that the slow time 



Schrodinger equation in Eq. (Al) has the identical form 



as the usual Schrodinger equation, and therefore, many 
results for static systems can be directly extended to pe- 
riodically driven syste ms t hrough the extension of the the 
inner product to Eq. (A4). 



2. Floquet Fermi golden rule 



" Floquet Fermi golden rule" gives the intuition behind 
the response of periodically driven systems under the in- 
fluence of perturbations. In particular, the result shows 
that quasi-energy of the effective Hamiltonian i^eff is a 
conserved quantity up to the driving frequency Q. In the 
context of our paper, this result implies electrons in the 
lower band of i^eff cannot be excited to the upper band if 
the frequencies of the perturbations, such as phonons or 
probe light, are smaller than the gap A. Thus the trans- 
port property of a non-equilibrium system described by 
i^eff is robust against phonon interactions as long as the 
chemical potentials of leads lie in the gap and phonon en- 
ergies are smaller than the gap A. Moreover, if one can 
prepare the system in the state with filled lower band of 
i^eff, then the gap of i^eff can be probed by observing the 
transmissions of low frequency lasers. 

This "Floquet Fermi golden rule" can be easily ob- 
tained through the "two-time" formalism described in 
the previous section. In the following, we consider the 
perturbations of the system such as phonons with fre- 
quency uj much smaller than Q such that uj <^ Q. We 
take the perturbation in the form V{r) = Ve~'^^^ . The 
usual derivation of Fermi golden rule can be applied in 
a straightforward fashion, and we obtain "Floquet Fermi 
golden rule", which gives the rate 7i^/ of exciting the 
initial Floquet state \(pi) to the final Floquet state \(pf) 
in the presence of the perturbation V{r); 



E mT\V\^'^))\'m -^uj-Ef- run). (A5) 



Here Ei and Ef are the quasi-energies of the initial and 
final Floquet states, respectively. In order to derive the 
result above, we represented the Floquet state \(pi) by the 
specific periodic state \^^). This choice is arbitrary and 
any other choice gives the same result. Since the physical 
Floquet state |(pf) can be represented as the states l^J^) 
for any integers m, the total transition rate is given by 
the sum of the rate from the state |<I>^) to states l^J^). 

This rate has the same form as the conventional Fermi 
golden rule, except for the summation over the Floquet 
energy index m. The delta function in the equation above 
imposes the conservation of quasi-energy which is the 
eigenenergy of effective Hamiltonian i^eff, which means 
the energy is conserved up to the driving frequency Q. 
This is a natural consequence of the fact that the system 
can absorb or emit the energy Q from the periodic drives. 

From this result, it is clear that such conservation of 
quasi-energy prevents the excitations of electrons from 
lower band to upper band when phonon energy u is 
smaller than the gap of the system, and Q is much larger 
than the total band-width of electrons. 
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3. Floquet adiabatic theorem 

In this subsection, we show, in analogy with the adi- 
abatic theorem of static systems, that a Floquet state 
follows an adiabatic change of Hamiltonian and stays in 
the Floquet state of the instantaneous Hamiltonian. This 
result indicates that the adiabatic increase of the inten- 
sity of light can be used to prepare the state with filled 
lower band of i^eff, whose properties can then be probed 
through low frequency lasers as argued above. 

St artin g from the slow time Schrodinger equation in 
Eq. (Al), we can follow the derivation of adiabatic the- 



orem and prove the analogous theorem for periodically 
driven systems. Here we briefiy outline the derivation. 

Suppose that the total Hamiltonian H(r) is slowly 
varying as a function of r. We are interested in how 
a Floquet state of H{0) at time r = evolves under this 
time evolution. Let \g) be the initial Floquet state and 



|G(ro)) be the result of evolving \g) under H{r) for time 
To- 

We denote the instantaneous eigenstates of H(r) as 
|a(r)) such hat H(r)|a(r)) = Ea{r)\a{r)). Then we 
express the state \G{r)) in terms of |q^(t)) as 



-i [ Eg{r')dt' 
Jo 



|G(r)) = exp -i 



c.(r)|^(r)) + ^c.(r)Kr))| (A6) 

In the absence of degenerate states, we can solve for 
the coefficients Cair) in the lowest order for the slow 
change of H amiltonian H(r) in the Schrodinger equation 
of Eq. (Al). The result is given by 



\G{r)) 



-i!^ Eg{T')dt' 



exp 



((«(T)||:|g(r))) 
Ep{T)-Eg{T) 



(A7) 



Thus to the zeroth order for the slow change of Hamil- 
tonian H(r), \G{t)) is the Floquet state of the instanta- 
neous Hamiltonian H(r) with possible accumulations of 
dynamical and Berry phases. The first order correction 
is given by the second term of Eq. (A7). 



For static systems of Dirac Fermions studied in this pa- 
per, we have shown that a gap proportional to opens 
at the Dirac point upon the application of light. If the 
chemical potential lies at the Dirac point before the ap- 
plication of light, the result above implies that the adia- 
batic increase of the intensity of light A{t) can be used to 
prepare the system close to the filled lower band state of 
i/eff- At exactly the Dirac points where the spectrum be- 
comes degenerate, the adiabatic theorem above does not 



apply, but these points represent only a tiny portion of 
the total states, and thus can be ignored for the calcula- 
tions of physical quantities. When the initial system is at 
finite temperature, such adiabatic increase of A{t) leads 
to non-thermal distributions of electrons in the spectrum 
of i^eff, but nonetheless the resultin g de nsity matrix can 
be calculated through the result Eq. ( A7 ) in the adiabatic 
limit. 

This Floquet adiabatic theorem can be used to obtain 
the Kubo's formulgP^in the non-interacting, periodically 
driven systems. In Ref.[9^, such result is applied to derive 
the extension of TKNN formula^ to periodically driven 
systems in infinite systems. 
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